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Abstract 

For the multisensor autoregressive moving average (ARMA) 
signals with unknown model parameters and noise 
variances, using the recursive instrumental variable (RIV) 
algorithm, the correlation method and the Gevers-Wouters 
algorithm with a dead band, the fused estimators of 
unknown model parameters and noise variances can be 
obtained. Then by substituting them into optimal fusion 
signal smoother weighted by scalars, a self-tuning 
distributed fusion Kalman smoother is presented. Using the 
dynamic error system analysis (DESA) method, it is 
rigorously proved that the self-tuning fused Kalman signal 
smoother converges to the optimal fused Kalman signal 
smoother, so that it has asymptotic optimality. A simulation 
example shows its effectiveness. 
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Introduction 

Information fusion techniques combine data from 
multi-sensors and related information to obtain more 
specific inferences than those obtained by using a 
single sensor [Liggins et al. (2009)]. The existing 
information fusion Kalman filtering is mainly focused 
on the information fusion Kalman filtering with 
known model parameter and noise statistics. However, 
in many applications, the model parameters and/or 
noise variances are usually unknown. The filtering for 
the systems with unknown model parameters and/or 
noise variances is called self-tuning filtering [Deng 
(2012)]. 

Recently, several self-tuning weighted fusion filters 
[Deng et al. (2008), Deng & Li (2007), Sun (2007)] have 
been presented, whose drawbacks are that only the 
noise variances are assumed to be unknown, while the 
model parameters are assumed to be known. In the 
existing results, most of them are for the single- 
channel signal systems with white measurement noise 
[Tao et al. (2010), Liu & Deng(2010), Gao & Deng 



(2010)]. And, up to the present, the information fusion 
is mainly focused on filtering fusion for the single- 
channel multisensor systems with white measurement 
noises, while smoothers for ARMA signals are rarely 
reported[Tao & Deng (2010), Tao & Deng (2012), Deng 
&Li (2007)]. 

In this paper, for multisensor autoregressive moving 
average (ARMA) signals with a common disturbance 
noise, using the Kalman filtering method, a self-tuning 
distributed information fusion Kalman smoother 
weighted by scalars is presented for the case when 
model parameters and noise variances are both 
unknown. By the dynamic error system analysis 
(DESA) method [Ran et al. (2009)], it is proved that the 
self-tuning distributed fusion Kalman smoother 
converges to the optimal distributed fusion Kalman 
smoother with probability one, so it has asymptotic 
optimality. Compared with[Tao & Deng (2010), Tao & 
Deng (2012), Deng & Li (2007)], the paper uses 
different signals model with a common disturbance 
noise. In addition, unknown model paremeters and 
noise variances are different from [Tao & Deng (2010), 
Tao & Deng (2012), Deng & Li (2007)]. 

Problem Formulation 

Consider the multisensor single -channel ARMA signal 
with L -sensor 

A{q- V )s(t) = C(q- l )w(t) (1) 
y,{t) = s{t) + v,(t) ^ i=l,2,--;L 

v,Q) = £(!) + €,(!) (3) 
where t is the discrete time, s(t) is the ARMA signal to 
be estimated, y t (t) is the measurements of the ;'th sensor, 
w(t) and v ; (f) are the input noise and measurement 
noises of the;'th subsystem, respectively, v,.(f) contain a 
common disturbance noise %(t) , q~ l is the backward 
shift operator, q~ l s{t) = s(t-l) , A(q~ h ) and C(<? -1 ) are 
stable polynomials having the form as A(q~ 1 )=l + 
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atf^ +---a„ a q~" a and C{q~ l ) = 1+ c x q 1 + ■■■c„ c q " c , n a is 
the order of polyno-mial A{q~ l ) , n c is the order of 
polynomial of C(q~ l ) . 

Assumption 1 w(t) , %{t) and e t {t) (z =1,2,---,L) are 
uncorrelated white noises with zero mean and 
variances o- 2 ,a 2 and a 2 , respectively, and the input 
noise w(t) and the measurement noise v,.(f) have the 
relation 



w(t) 
vXt) 



r (0 vj(0] 



< 



(4) 



where the symbol E denotes the mathematical 
expectation, the superscript T denotes the transpose, 
S n =l,S tk =0(t*k). 

From (3), it can be seen that the measurement noises 
v t (t)(i = 1,- -,L) are correlated,i.e. 

R tj = E[ v, (0 v ] {t)] = a] (i *j),R,= E[v, (t)vj (f )] = a 2 + a 2 (5) 

Assumption! A(q x )snd C(g _1 )are stable polynomials. 

Assumption 3 The model parameters fl , , c , , and noise 
variances <j 2 w ,<j\ and a t J are unknown. 

Assumption 4 The realizations of measurement 
stochastic process y t (t) (i = 1,2,---,L) are bounded for f, 
with probability one. 

The Optimal Fusion Kalman Signal Smoother 
Weighted by Scalars 

Setting w(t ) = w(f - 1) , yields that w(t ) has the variance 

The ARMA signal (1) can be rewritten as 

A{q- l )s(t) = C(q- l )w(t) (6) 
C (q' 1 ) = c iq ~ [ +■■■ + c n _q~ n " , c i = c w , j = 1,2, • • • , 

M C -=M C +1 (7) 

From (2) and (5), we have the state space model with 
the companion form 

x(t + 1) = 0x(t ) + Tw{t) (§) 

y,(f) = #x(f) + v,(f), i = l,2,-,L (9) 

s(t) = Hx(t) (10) 

-a, 



"o-l 

-a o ••• o 



//=[l o ■■■ o] 

r T = [l Cl - cJ,Cj=<Kj>n e ) 



Lemma 1 [Deng (2012)]. For the multisensor system 
(7)-(9) with known model parameters and noise 
variances, the /th sensor subsystem has the local 
steady-state optimal Kalman predictor 

x i {t + \\t) = >F Pi {t)x i {t\t-\) + K Pi {t)y i {t) 
r pi (t) = 0-K Pi (t)H 

K Pi (t) = 0Z i (t\t-V)H T Q- 1 (t) 



(12) 
(13) 

(14) 

QJt) = Hi;,(t\t-l)H T +R, ( i5) 

where the local prediction error variance matrix 
Z t {t 1 1 - 1) satisfies the Riccati equation 

27 (t + 1 1 1 ) = ,{t 1 1 - 1) - 27 (f 1 1 - l)H 1 (HS l (t 1 1 - 1) x 
H 1 +R l ) l HZ,{t\t- \)]0 T + ralr 1 (i 6) 

and the local prediction cross-covariance 

27 (t 1 1 - 1) = E[3c ; (f If- l)5c / (f I f - 1)] , 

with x t (t 1 1 - 1) = x(f ) - X; (? I f - 1) satisfies the Lyapuno v 

equation 

27 (f + 1 1 1 ) = !P„ (f)27. (t\t- (r) + /^T T + 

K pi (t)RijKpj(t) i*j i,j = \,2,-,L (17) 

Lemma 2 [Deng (2012)]. For the multisensor system (1) 
and (2) with known model parameters and noise 
variances, the ! ™ sensor subsystem has the local 
optimal Kalman signal smoother as 

JV 

x.(t-N\t) = x i (t-N\t-N-\) + 'Y d K i (t-N\t-N+ j)x 

e,(f-N+j),i = l — ,L (is) 
The optimal Kalman signal smoother is given by 

s i (t-N\t) = Hx i (t-N\t) (19) 
where the time-varying optimal Kalman predictor 
x t (t- N \t- N-l) can be computed by (12), and we 
have 

s,(t) = y l (t)-Hx l (t\t-\) (20) 

7-1 

K t (t\t + j) = Z t (f I f - l){fj Tj,(t + k) }H T Q,' (t + j) (21) 

and the cross-covariance among smoother errors 
l i (t-N\t) = s(t-N)-s i (t-N\t) / i s P sij (t-N\t) = HP ij (t-N\t)H T , 

and the error variance matrices and covariance 
matrices of local Kalman smoother are given as 

N 

P i (t-N\t) = £ i (t-N\t-N-l)-'J]K.(t-N\t-N + j)x 

QAt-N + j)Kj(t-N\t-N + j) (22) 
PM-N\t) = r„(!-N)ZM-N\t-N-Wl(!-N) + 



(11) 



Ifc< 

p=0 



^< p (t-N) KJ p (t-N) 





'Kf p {t- 


N) 


_0 R tJ _ 




N)_ 



(23) 
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with the definition P.. (t - N\ t) = P, (t - N 1 1) . 

N 

<F lN (t-N) = I„-Y J K l (t-N\t-N + k)Hx 
>F pi (t-N + k,t-N) . 

N 

K^{t-N)=-^K i {t-N\t-N + k)Hx 

k=p+\ 

*¥ pi {t-N + k,t-N + p + \)r p = 0,---,N-V, 

f 

N 

K v ip {t-N)= ^K.(f-N\t-N + k)Hx 

k=p+l 

r pi (t-N + k,t-N + p + l)K pi (t -N + p)- 
K.(t-N\t-N + p) p = 0,--,N-l; 

(t-N) = O t KJ N (t-N) = -K, (t-N\t) (24) 
where ¥ pi (t + k,t) = i F pi (t + k-\)---'F pi (t), W pi (t,t) = I„. 

Lemma3 [Deng (2012)]. For the multisensor system (7)- 
(9) with the assumptions 1 and 2, the optimal 
information fusion Kalman signal smoother 
s (t - N \t) is given by 

L 

s (t -N\t) = ^eo i (t-N\ t)s t (t -N\t) ^ 

where the optimal scalars weighting coefficient vectors 
a>i(t-N\ t) is given by 

[a\(t-N\t), a) L (t-N\t)]=[e T (P; 1 (t-N\t)f [ x 



er l e J (P(t-N\t)y 



(26) 



where e T = [1, • • • ,1] , defining the LxL matrix as 

P s (t-N\t) = (tiP sij (t-N\t)), i,j = l,2,-,L (27) 

Self-tuning Fusion Kalman Smoother Weighted 
by Scalars 

When model parameters and noise variances are 
unknown, substituting estimators into the optimal 
fusion Kalman smoother will yield the self-tuning 
fusion Kalman smoother which consists of the 
following steps: 

Step 1. Applying the system identification algorithm 
[Ljung (1999)], the estimator a t (t) can be obtained. 
Applying the correlated method [Ran et al. (2009)] and 
Gevers-Wouters algorithm [Gevers & Wouters (1978)], 



R t {t) can be obtained, where Ry(t) = ct^O) (i * j) , 

R f (t) = <t|(0 + 0^(0 . And the model parameters and 
noise variance estimators are consistent, i.e. 

a, (t) -> a t , c, (t) -» c, a 2 w (r) -> a 2 w , a 2 (t) -> a 2 , 



a\(t)^>a\ R i (t)^R i Ry(t)^Ry 
as t -> oo, w.p.l (28) 
Step 2. Substituting the estimators a, (f) and c,(f) into 
(11) yields S(t) and f(t). In Lemma 1,0 r a 2 w R s 

/ / / 

and Rij {i, j = !,■■■, L) are replaced by <P(t) f(t) al(t) 

/ / / 

Rj(t) and Rq(t) , respectively. And the estimators 
E t (t \t-l) satisfy the self- tuning Riccati equations 

E i (t + 1 1 1) = 0{t)[Z i (t 1 1 - 1) - I i {t 1 1 - 1)H T (H x 

i i (t\t-\)H J + R t (?))"' HZ i (t\t- 1)]0 T (0 + 

f(t)a 2 w (t)f T (t) (29) 

and the prediction cross-covariance matrices satisfy 
the self-tuning Lyapunov equations 

i, (t + 1 1 1 ) = !P H (f )i, (f I f - 1)^ (0 + f(t)al(t)f T (t) + 

K pi {t)R tj (t)K T PJ (t) j, j = 1,2, • • • ,L (30) 

with the definition Z H (t + 1 1 1) = Z i (t + 1 1 1) . Substituting 
the estimators into (12), the self-tuning local Kalman 
predictor can be given as 

X* (t + 1 1 = W n (f )*,' (t \t - 1) + K Pi (t) y t (t) (31) 

V Pi (t) = 0(t)-k pi (t)H (32 ) 

k Pi (t) = S(t)i i (t\t-l)H T Q- 1 (t) (33 ) 

Q^t) = H£,(f\t-l)H ,I +R l (f) (34) 

Substituting the estimators into (18)-(24), yields the 
self-tuning local Kalman smoother x- (t-N\ t), 
(i = 1, • • • , L) as 

x,:'(f-^lf) = -J/(r-A ? lf-A ? -l) + ^^,(r-A ? lr-A ? + j)x 



e i (t-N + j),i = L 



ff f (r) = ;y,(f) I f-1) 
f ,(f I f + J) = Mt 1 1 -l){Yl^(f + k)}H T (H x 



(35) 
(36) 



£,(.t + j\t+j-l)H T +R,(t)r (37) 

The error variance matrices and covariance matrices of 
self-tuning local optimal Kalman smoother are given 



fusion estimators c,(f), a 2 w (t) cr f (0 a 2 (t) Ry(t) and as 



P i (t-N\t) = i i (t-N \t-N-l)-^ j k i (t-N\t-N+ j)x 



7=0 



Q d (t-N + j)Kj(t-N\t-N + j) 



(38) 



^.(f-/vif) = ^(f-^(f-^vif-^v-i)^a-^v) + 
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f\K;(t-N) Kj p (t-N) 

p=0 



Kit) o 

RM) 



KjJ(t-N) 



(39) 



(40) 



with the definition P u (r - N I f) = P, (f - N I f) . 

(V 

'F iN (t-N)=I n -^K i (t-N\t-N + k)Hx 

>F pi {t-N + k,t-N) . 
£ £ ( f - A?) = - ^ £ ,( f - ^ 1 1 - N + k)H x 

k=p+l 

<f pi (t-N + k,t-N + p + l)f(t) p = 0,-,N-ii 

/ 

N 

Kj p (t-N)= y* i K t (t-N\t-N + k)H x 

<F pi {t-N + k,t-N + p + \)K pi {t-N + p)- 

K t (t-N\t-N + p) p = 0,---,N-l; 

K£(f - N) = 0K* N (t - N) = -K,(f -N\t) 

where *P pi {t + k,t) = 9 r pi (t + k -!)■■■ <P pi (f ), <P pi (t, t) = I n . 

And the self-tuning local Kalman signal smoother is 
given by 

s° (t - N 1 1) = Hx< (t-N\t) (41) 

Step 3. Applying (26) and (27), the estimates 
P s (t-N\t) and & i (t-N\t) are obtained, and the self- 
tuning fused Kalman smoother is given as 

L 

S; (t-N\t) = ^d> i (t-N\t)s°(t-N\ t) ( 42 ) 
The above three steps are repeated at each time t. 

Convergence Analysis 

Convergence of Self-tuning Lyapunov equation 
Theorem 1. For the multisensor system (1) and (2) with 

the assumptions 1-4, the solution ^(f + HO of the 
self-tuning Lyapunov equation (30) converges to the 
solution 2q(t + 11?) of the optimal Riccati equation (17), 
i.e. 

[4(f + llf)-^(f + H0]^0,as f ^oo,i.a.r (43) 

Proof. It has been proved that the solution Z^t + Wt) 
of the self-tuning Riccati equation converges to the 
solution Z i of the steady-state optimal Riccati 
equation [Tao & Deng (2010)], i.e. 

Z i (t + l\t)-Z i ;-»0,as f^oo, i.a.r (44) 

According to the stability theory of the Kalman 
filtering [Kamen & Su (1999)], the solution Z.(t + l\t) 
of the optimal Riccati equation converges to the 



solution Z ( of the steady-state optimal Riccati 
equation , i.e. 

z i {t+\\t)-i i ->o (45) 

From (44) and (45), we have that 

Z i (t + l\t)-Z i (t + l\t) = {Z i (t + l\i)-Z i )- 
(i7 i a + H0-^ i )^0,as f->oo, i.a.r (46) 

From (14) and (33), according to the stability theory of 
the Kalman filtering [Tao & Deng (2010)], the steady- 
state optimal Kalman predictor gain matrix K Fj i s a 
continuous function with respect to , Z j and R i , i.e. 
K Pi = 0Z i H T (HZ i H T +RJ 1 =f l (0,Z l ,R i ) . From (11) 
and (28), we have &(t)->0,f(t)->r. According to 
the continuity of f , applying (28) and (46) yields that 

[K Pi (t) - K PI (f)] -> , as t -> oo , i.a.r (47) 

Applying (13), the steady-state predictor transition 
matrix is a continuous function with respect to K pi 
and<£, i.e. !P„ = <P-K Pi H = g,{K n ,^) . From (47) and 
the continuity of g t , we can yield that 

\W„(t) - !P„(0] -> 0, as t -> oo , i.a.r (48) 

Setting &„(t) = W pi (t) + A &„(t) , K Pi (t) = K Pi (t) + AK pi (t) , 
then 

A W Pi (t) -> 0, AK Pi (t) -> , as t -> oo , i.a.r (49) 

Subtracting (17) from (30), and setting 

E ij (t)=Z ij (t + l\t)-Z ij (t + l\t) (50) 

we obtain the dynamic variance error Lyapunov 
equation 

Ey(t) = Y^E^t- l)iF£ (*)+£/„(*) (51) 

with the input 

U t (0 = A <F Pi (f)£ 9 (t\t- 1)5P* (t) + A W Pi (r)i, (f 1 1 - 1) x 

A 4>l it) + <F Pl (t)£ 9 (t\t- 1) A $1 (t) + f(t)al(t)f T (0 - (52) 

ra 2 w r T + K Pi {t)R tj (t)K T P] (f) - K Pi WyKjj (0 

According to the stability theory of the Kalman filte- 
ring [Kamen & Su (1999)], applying Lemma 6 yields 
that K Pi {t), Ypiit), fiptit) and K Pi (t) are uniformly 
asymptotically stable. Applying the boundedness of 
f(t) , &l(t) yields that f(t)a 2 w (t)f J (t) is bounded. 

Similarly, it can be obtained that K Pl (t)Rjj(t)Kj,j(t) is 
bounded. Applying Lemma 3 [Tao & Deng (2012)] to 
(30) yields Z tj (t 1 1 -1) is bounded. Therefore, from (52), 
applying (28), (43) and (49) f(t)^r yields that 
Uyit) -> 0, as t -> oo , i.a.r . Applying Lemma3 [Tao & 
Deng (2012)] to (51) yields 
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£, y (0^0,as f->oo,i.a.r (53) 
which gives (43). The proof is completed. 

Convergence of the Self-tuning Local And Fused Kalman 
Smoothers 

Theorem 2. For the multisensor systems (1) and (2) 
with the assumptions 1-4, the local self-tuning Kalman 
predictor x-(t + l\t) converges to the local time-varying 
optimal Kalman predictor x ( (f + 1 If) in the sense that 

[x- (f + 1 1 f) - x t (f + 1 1 f)] -> 0, as 1 -> 00 , i.a.r. (54) 

Proof. Denoting J. (f ) = x. (t + 1 1 1) - x, (f + 1 1 f) , subtracting 
(12) from (31) yields a dynamic error system 

5 i (f) = 9 r „(f)S t (f-l) + u t (f) (55) 

u, (f) = zl !F„ (Ox/ (f I f - 1) + ^ (f )y ( (0 (56) 

From the boundedness of K Pi (t) and the assumptions 
4, K Pi {t)y t {t) is bounded. Since ¥p,-(f) is uniformly 
asymptotically stable, hence, the application of Lemma 
4 [Tao & Deng (2012)] to (31) yields that x/(f I f -1) is 
bounded. From the boundedness of x* (f I f - 1) and 
y t {t) applying (49) to (56) yields that 

«;(f) -> , as f -> oo , i.a.r. (57) 
Since f Pi (f) is uniformly asymptotically stable, 
applying Lemma 4 [Tao & Deng (2012)] to (55) yields 
that 

Sf(t) -> , as t -> oo , i.a.r. (58) 
Hence, (54) holds. The proof is completed. 

Theorem 3. For the multisensor system (1) and (2) with 
the assumptions 1-4, the self-tuning fused Kalman 
smoother x s Q (t-N If) converges to the optimal fused 
Kalman smoother x (t-N\t) m a realization, i.e. 

[jco(f-AMf)-jc (f-Anf)]->0, as f->oo, i.a.r. (59) 

Proof. From (20) and (36), applying (54) yields 

[3(f) -£,.(*)]-»(), as f-> oo, i.a.r (60) 

From (21) and (37), applying (28), (43) and (48) and the 
property of the continuous function yields 

[K i (t-N\t-N + j)-K,(t-N\t-N + ;')]-)-0 7 

as t oo , i.a.r (61) 
Setting s. (f ) = s t (f ) + Ae i (t) ,K i {t-N\t-N + j) = 
K i (t-N\t-N + j) + AK i (t-N\t-N + j), from (60) and 
(61), we have 

As,(t) — »0, AK i (t — N\t — N + j)—>0 

as t oo , i.a.r (62) 
Denoting d t (t) = x- (f - N I f) - x t (f - N I f) , and subtractting 
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(18) from (35) yields 

J l .(0 = jc, J (f-Af|f-A^-l)-x i (f-Af|f-Af-l) + 7 i (f) (63) 

JV 

rj,(t) = Y J K i (t-N\t-N + j)s t (t -N + j)- 

K^t-Nlt-N + j^it-N + j) (64) 
From (64) we have 

/7,(0 = Y,K t (t-N\ t-N + j)As i (t-N + j) + 

j=o 

AK i (t-N\t-N + j)e t {t -N + j) + 
AK i (t-N \ t-N + f)Ae t (t-N + j) (65) 
From (36), applying the boundedness of y,(f) and 
x s i {t\t-\) yields that £,(f) is bounded. From (65), 
applying (62) yields 

77,-(f)->0, as t -> oo , i.a.r (66) 
From (54) and (65), we have 

d ( (t) ->■ 0, as f -» oo , i.a.r (67) 

From (22), (23), (38),(39), applying (28), (46)-(48) and 
(61), we have 

P.{t-N\t)^P t (t-N\t) P^{t-N\t)-^P v {t-N\t) 

as t ->oo, i.a.r (68) 
From (26) and (27), each element of (o,(t-N\t) j s a 
continuous function of elements of P kj (t-N\t) / we 
have 

a> i (t-N\t)^co i (t-N\i), as f^oo, i.a.r (69) 

Setting co^t-N \t) = co i (t-N\t) + Aa> t (t-N\t), from (69) 
we have that zl<» ; (f -iV 1 1) -> , subtracting (25) from 
(42) yields 

Jcq (f - ATI f) - ic (f - N I f ) = ^ fi) ; (f - AT I f ) x 

L 

[r; (t-N\t)-x i (t-N\t)] + ^ J Aa> i (t-N\ t)x° (t-N\t) (70) 

Applying (69), zl<^ (t - N\ t ) -> , and boundedness of 
x\ (t-N\t) yields that (59) holds. 

Simulation Example 

Consider the ARMA signal s(t) with measure-ment 
noises 

A(q- 1 )s(t) = C(q- 1 )w(t) (71) 
y i (t) = s(t) + Z(t) + e i (t), i = \,2,-,L (72) 
where £(f) , e,(f) and w(f) are independent Gaussian 
white noises with zero mean and variances <r|, a 2 ei 
and cj\ . In simulation, we take 

A(q- ] ) = 1 + a^ 1 + a 2 q- 2 = 1 + Olq' 1 - 0.48g~ 2 , 
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C{q- i ) = l + c l q- 1 



= 0.1 ,o* =0.7, 



< = 0.1, cr e 2 2 =0.2, cr; 3 = 0.3. 
The problem is to design the self-tuning information 
fusion filter of s(t) , where cr 2 , cr 2 , cr 2 . (;' = 1,2,3) , a, , a 2 
and c, are unknown. 

The simulation results are shown in Figure 1-Figure 6. 
From Figure 1-Figure 5 in which the curved lines 
denote the estimates, and the straight lines denote the 
true values, it can be observed that the estimates of the 
model parameters and noise variances are consistent. 
Figure 6 shows that the self-tuning fusion Kalman 
signal smoother Sg(f-lU) converges to the optimal 
fusion Kalman signal smoother s (t -1 1 1) , i.e. its error 

converges to zero. 
i 
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FIG. 3 THE FUSED ESTIMATE CURVE OF NOISE VARIANCE 
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FIG. 4 THE FUSED ESTIMATE CURVES OF NOISE VARIANCES 
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FIG. 5 THE FUSED ESTIMATE CURVES OF MODEL PARAMETER 
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FIG. 6 THE ERROR CURVE BETWEEN THE SELF-TUNING AND 
OPTIMAL FUSION KALMAN SIGNAL SMOOTHERS 

e(t) = F (t-l\t)-s (t-l\t) 

Conclusion 

For multisensor ARMA signal with unknown model 
parameters and noise variances, by the system 
identification, correlation method and the Gevers- 
Wouters algorithm with a dead band, a self-tuning 
distributed information fusion ARMA signal Kalman 
smoother was presented. By the DESA method, it was 
rigorously proved that it converges to the optimal 
fusion ARMA signal Kalman smoother with 
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probability one, so that it has asymptotic optimality. 
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